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Abstract 

The bulk matter Randall-Sundrum (RS) model is a setup where Standard Model (SM) 
matter and gauge fields reside in the bulk of 5D warped spacetime while the Higgs field is 
confined on the IR brane. The wavefunctions of the 1st and 2nd generation matter particles 
are localized towards the UV brane and those of the 3rd generation towards the IR brane, so 
that the hierarchical structure of the Yukawa couplings arises geometrically without hierarchy 
in fundamental parameters. This paper discusses observing signals of this model in the case 
where the Kaluza-Klein scale is far above the collider scale, but the model is combined with 
5D Minimal SUSY Standard Model (MSSM) and SUSY particles are in the reach of collider 
experiments. A general SUSY breaking mass spectrum consistent with the bulk matter RS 
model is considered: SUSY breaking sector locates on the IR brane and its effects are mediated 
to 5D MSSM through a hybrid of gravity mediation, gaugino mediation and gauge mediation. 
This paper argues that it is possible to observe the signals of the bulk matter RS model through 
rare decays of "almost SU(2) singlet mass eigenstates" that are induced by fiavor- violating 
gravity mediation contributions to matter soft SUSY breaking terms. 



1 Introduction 



The origin of the fermion mass hierarchy is a long-standing mystery in particle physics. In the 
Standard Model (SM), it is explained by the hierarchical Yukawa coupling constants, but the 
hierarchy itself is still introduced by hand. Many authors have proposed models beyond the SM 
where the Yukawa coupling hierarchy arises from non-hierarchical couplings of a fundamental 
theory. The bulk matter Randall- Sundrum (RS) model [H [2] is one of the successful models. 
In this model, SM fermions are identified with the zero-modes of 5D Dirac fermions that live 
in 5D warped spacetime (bulk), whereas the Higgs field is confined on 4D spacetime of the 
infrared (IR) brane. With non-hierarchical values of 5D Dirac masses, the zero-modes of the 
5D fermions can be localized towards either the ultraviolet (UV) brane or the IR brane. In 
this way, the geometrical overlap between each zero-mode in the bulk and the Higgs field on 
the IR brane gains exponential hierarchy, which gives rise to the hierarchical structure of the 
Yukawa coupling constants. As far as the model gives a natural explanation to the fermion 
mass hierarchy, neglecting the gauge hierarchy problem, it is sufficient that the warp factor 
be around ~ me/ nit, or equivalently the Kaluza-Klein (KK) scale be around ~ M* ■ {me/mt) 
(M,, indicates the 4D reduced Planck mass) if the 5D Planck scale is the same as that of 4D. 
Then the KK modes need not exist at TeV scale. Futhermore, the physics of flavor has already 
imposed severe constraints on the mass of KK modes; for example, the 1st KK gluon in the 
bulk matter RS model induces flavor changing neutral current interactions. The data on the 
— mixing require that its mass be larger than 21 TeV [4]. 

If the KK modes appear only far above the TeV scale, it is impossible to produce them at 
colliders and conflrm the model. In this paper, I argue that it is possible to observe indirect 
signatures of the bulk matter RS model at near-future colliders. I consider the case where 
the KK scale is at an intermediate scale between Planck and TeV, but the 4D effective theory 
contains M = 1 SUSY which is broken at TeV scale and it can be described by Minimal SUSY 
Standard Model (MSSM). Hence the SUSY particles are accessible at colliders while the KK 
modes are beyond their reach. This is a natural situation because SUSY breaking at TeV 
scale is necessary to solve the gauge hierarchy problem when the KK scale is at an intermediate 
scale. I consider a general SUSY breaking mediation mechanism that is consistent with the bulk 
matter RS setup, in contrast to the paper [5] where a simultaneous explanation to the SUSY 
breaking mediation mechanism and the Yukawa coupling hierarchy based on the RS spacetime 
was pursued. In the setup of this paper, SUSY breaking sector locates on the IR brane and 
its effects are mediated to MSSM in the bulk through contact terms on the IR brane (gravity 
mediation [6]), renormalization group evolutions below the KK scale (gaugino mediation [7]) 
and gauge interactions with messenger flelds on the IR brane (gauge mediation [8]). (Anomaly 
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mediation contributions [9] are suppressed at least by the warp factor compared to gaugino 
mediation ones and hence are neghgible.) 

Gravity mediation contributions are the key to observe signals of the bulk matter RS model. 
This is because they arise from contact terms on the IR brane in a similar way to the Yukawa 
coupling constants. The basic strategy for testing the model is as follows. Since the zero-modes 
of matter hypermultiplets (SM fermions are their fermionic components) reside in the bulk, 
they have contact terms with SUSY breaking sector on the IR brane which induce soft SUSY 
breaking terms through gravity mediation. The amount of the gravity mediation contribution 
to each term is proportional to the geometrical overlap between the zero-mode fields and the 
IR brane. In the bulk matter RS model, the same overlap also gives rise to the hierarchy of 
the Yukawa coupling constants. Hence the flavor structure of gravity mediation contributions 
and the Yukawa coupling hierarchy are related and one can predict the former from the lat- 
ter. (A similar setup |T0] was proposed on a different context as a 5D realization of "flavorful 
supersymmtery" [H].) Therefore the model can be tested through the measurement of gravity- 
mediation-origined soft SUSY breaking terms. Our task is then to extract gravity mediation 
contributions from experimentally observable quantities related to soft SUSY breaking terms. 
To be a realistic SUSY breaking model, the model must contain dominant flavor-conserving soft 
SUSY breaking mass terms that arise through gaugino mediation and gauge mediation. On the 
other hand, gravity mediation intrinsically violates flavor. Therefore one can indirectly measure 
the gravity-mediation-origined terms through flavor- violating interactions of SUSY matter par- 
ticles. One obstacle is that the Yukawa coupling constants themselves induce flavor-violating 
soft SUSY breaking terms through renormalization group (RG) evolutions, as in models with 
minimal flavor violation. However, one can distinguish gravity mediation contributions from 
the RG effects of the Yukawa couplings by focusing on SU(2) singlet SUSY matter particles, for 
which flavor violation of gravity mediation contributions can be signiflcantly larger than that 
of the RG effects. In this paper, I introduce three promising channels for observing signatures 
of the bulk matter RS model at colliders. One is that SU(2) singlet smuon mixes with stau 
through gravity-mediation-origined soft terms and one measures the branching ratio of "almost 
singlet smuon mass eigenstate" decaying into SM tau and another SUSY particle. Another 
channel is that SU(2) singlet smuon mixes with selectron and one measures the branching ratio 
of "almost singlet smuon mass eigenstate" decaying into SM electron and another SUSY par- 
ticle. The third is that SU(2) singlet scharm mixies with stop and one measures the branching 
ratio of "almost singlet scharm mass eigenstate" decaying into SM top and another SUSY par- 
ticle. The bulk matter RS model predicts these branching ratios and can be tested through 
their measurements. In this setup, one cannot determine the exact values of the contact term 
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couplings of the 5D theory. I here assume that the contact term couphngs for gravity media- 
tion are all 0(1), and estimate the orders of magnitudes of flavor- violating soft SUSY breaking 
terms. Still the bulk matter RS model gives predictions on the magnitudes of the branching 
ratios of sparticle rare decays. 

This paper is organized as follows. In section 2, I review the bulk matter RS model and 
combine it with 5D MSSM. 5D disposition of matter wavefunctions are determined so as to 
reproduce the hierarchical structure of the Yukawa coupling constants, Cabbibo-Kobayashi- 
Maskawa (CKM) matrix and the neutrino mass matrix with 0(1) couplings of the 5D theory. 
In section 3, I derive the SUSY particle mass spectrum with emphasis on its flavor-violating 
sector. In section 4, I make predictions on the flavor mixings of SUSY particles, and compare 
the bulk matter RS model with other models. In section 5, I discuss experimental methods to 
test these predictions. 



2 The Bulk Matter RS Model with 5D MSSM 

2.1 Setup 

Consider 5D warped spacetime with the metric [1] : 

ds^ = e-^'^^y^ri^^dx^dx" - dy" , (1) 

where y is the 5th dimension compactified on the orbifold /Z2 ■ —irR < y < ttR , and k is 
the Anti-deSitter curvature that is of the same order as the 5D Planck scale M5. Assuming 
that the warp factor, e"^'^'^, is much smaller than 1, we have the following relation for k and 
Ms : 

Ml = i|^(l-e-^'=^-) ^ ^, (2) 

where is the 4D reduced Planck mass, which implies k ~ M5 ~ M^. The UV brane is 
put aX y = and the IR brane ai y = ttR. The Planck scale on the UV brane is M5, while 
that on the IR brane is M^e~^^'^ . The IR scale, ke~^^'^ ~ M^e~^^'" , is assumed to be at an 
intermediate scale between and TeV. In particular, we assume that it is far above 21 TeV. 
Since the most severe constraint on the IR scale of the bulk matter RS model comes from the 
data on the — mixing, which require it to be larger than 21 TeV my model is free 
from any constraint on the bulk matter RS model itself. At the same time, it is hopeless to 
observe the effects of the Kaluza-Klein excitations by near-future experiments. 

Consider 5D MSSM [3] where the 4D M = 1 Higgs superfields are confined on the IR 
brane, and the 5D M = 1 gauge superfields and matter hypermultiplets live in the bulk. In 
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the following, we use the 4D superfield formalism extended with the 5th dimension y. We 
introduce a chiral superfield, X, on the IR brane whose F-component, Fx-, develops vev to 
break 4D = 1 SUSY there. We consider both cases where there are one to several messenger 
pairs on the IR brane and there is no messenger pair at all. (It is easy to extend the model to 
cases where the messengers live in the bulk.) The gauge symmteries of the messenger pairs are 
not specified. SU(2) doublet squark, singlet up- type squark, singlet down-type squark, doublet 
slepton, singlet charged slepton hypermultiplets are denoted by Qi,Ui, Di, L^, Ei, respectively, 
with i being flavor index. The up-type Higgs and the down-type Higgs superfields are denoted 
by Hu, Hd, respectively. 

An off-shell 5D A/" = 1 gauge superfield consists of a 5D gauge field Am {M = 0, 1, 2, 3, 5), 
two 4D Weyl spinors Ai,A2, a real scalar S, a real auxiliary field D and a complex auxiliary 
field F, all of which transform as the adjoint reprsentation of some gauge group. They combine 
to form one 4D A/" = 1 gauge superfield V and one 4D A/" = 1 chiral superfield x that are 

X = -^{^ + iA5) + V29X2 + e^F . 
v2 

By Z2 : y —y symmetry, they transform as 

V ^ V , X ^ -X- 

Assuming the invariance of the theory under the Z2 symmtery, we obtain the following action 
for 5D A/" = 1 gauge superfields: 



S.Daauge = Uy d^X C'^'^l^ 



4(^?5") 



^ [ dVI^I tr| (ei'=l^liy"")(ei*^l^lH/") + h.c. 



(3) 



where a labels gauge groups and W"""^ denotes the field strength of V"' in 4D fiat spacetime. 
When the unitary gauge, = 0, is chosen, only V"" has a massless mode in 4D picture. This 
mode has no dependence on y and will be written as Vq{x, 6, 9). 

A 5D M = 1 hypermultiplet is expressed in terms of two 4D M = 1 chiral superfields $, 
that are in conjugate representaions of some gauge group. We assume that the former is Z2- 
even and the latter Z2-odd. Taking the basis of diagonal bulk masses, we have the following 
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action for 5D JV = 1 hypermultiplets: 



Si 



5D chiral 



J d^^e'^l^l $^{9y-x/V2- (3/2-Ci)A;}$i + h.c. 



(4) 



where i is a flavor index and q denotes the 5D bulk mass in unit of AdS curvature k. Only $j 
has a massless mode in 4D picture, which will be written as 6')e^^/^~'^'-*'^'2''. 

We write down the 4D effective action for the fields in the bulk in terms of the massless 
modes: 



SiDeff. — J ^ ^ 

+ fd^e2 



■ 27tR 

(1 - 2Ci)k 



d^e w^^w^ + h.c. 



(5) 



where the dimensionful 5D gauge coupling, g^, is connected to 4D gauge coupling by the 
relation: g^ = V^nRg^. 4>i represents the zero-mode of each of Qi, Ui, Di, Li, Ei. 

We also introduce an IR-brane-localized action. Below are the parts of the action relevant 
to the topic of this paper. 

MSSM term: 



Sir D / d^x 



I d'x I d'e e-''^- { Hie-'^H^ + H^H, } 



{yd) 



kl 



HdDkQi 



+ I d'9 e-''^""^ e(3-c^-c„)fcij. yy^HdE^Ln + h.c. 

Ms 



Gaugino mass term: 



Sir D I d X 



Matter soft SUSY breaking mass term: 

Sir D j d^x 

+ {Q ^ U, D, L, E) . 
A-term-generating term: 

Sir D / d^x 



d^e da^w^w;^ + h.c. 

Ms 



h.c. 

(6) 

(7) 



Ml 



QlQj 



Mi Mi 



Mi 



h.c. 



(9) 



Messenger term: 



+ I (^'9 e-=^'=^- { M^,,si ^i^i + Xmessi X EjEr } + h.c. 



(10) 



where Mmessi indicates the SUSY conserving mass for the messenger pair S/, S/. Note that we 
did not necessarily assume the existence of messengers. In that case, only gaugino mediation 
gives rise to flavor-conserving soft masses, as is realized in the model [5]. 

In addition, the terms for the Higgs superfields exist on the IR brane. We simply assume 
that yU-term and i?/i-term are somehow derived at TeV scale. 

We normalize X, if^, Qj, f/j, -Dj, Lj, Ei, Ej, to make their kinetic terms in the 
4D effective theory canonical. This is done by the following rescaling: 

X ^ X = e-^^'^X, Hu ^ Hu = e-'^^'^Hu, ^ = e-'^^'^Hd, 



^{l-2ci)kRTT _ I 

' {l-2ci)k 



— e — — e ^j. 



Then the MSSM term, the gaugino mass term, the matter soft SUSY breaking mass term, 
the A-term-generating term and the messenger term become as follows: 



Sir D j d^x 



j d^x j d'e [ Hle-^'H^ + HyH, } 



1 - 2c, 



1 - 2c, 



2{1 - e-(i-2c.)M^} y 2{1 - e-(i-2c,)fcR^} 

k 



+ 




2{1 - e-(l-2c^)fcR7r| V 2{l - e-{l-2c„)fcR7r| 

The gaugino mass term will be 



{yd)ki HdDkQi 



{ye)mn HdEm.Ln \ + h.C. 



.(12) 



Sir D / d^x 



d?e d. 



X 



(13) 
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The matter soft mass term will be 



Sir D / d^x 



1 - 2c, 



1 - 2c.- 



2{1 _ e-(i-2'=')'=^'^} y 2{1 - e-(i-2c^)'=«^} M5 



X 



+ ( g ^ f/, L, E ) . 



The A-term-generating term will be 



(14) 



Sir D 



+ 




XHJJiQj 



2{1 - e-(i-2cm)fcfl7r} y 2{1 - e-(i-2c")'=^'^} M5 Mgc-'^^^ 



The messenger term will be 



h.c. 



(15) 



d^'e { M^essie-^^^ EiEi + XmessiX EiEi } + h.c. 



(16) 



We introduce light neutrino masses by writing an IR-scale-suppressed higher dimensional 
operators or by adopting the seesaw mechanism [12]. In either case, we have the following term 
for light neutrino masses: 



Sir D y d'x y d^e e-'"''^ ^{3~c,^c,)kR.n ^y^-^^^ 



LpHuLqHu 
— M 



d^x / d^e 



1 - 2c., 



1 -2c„ 



2{1 — e-(l~2cp)fc_R.7r| y 2|1 _ ^-{l-2cq)kR^x^^ 



LpHuLqHu 

M 



kR-K 



h.c. 



(17) 



where Mgeesaw indicates the mass scale relevant to the light neutrino mass. To maintain the 
generality of the model, we hereafter consider cases with SU(2) singlet neutrinos. Their Majo- 
rana masses in the 4D effective theory are assumed to be around a common scale, denoted by 
^Maji which is lower than Mmessi^''^^^ or M^e~^^'" . The results of this paper can be extended 
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to cases without singlet neutrinos by dropping terms containing MmuJ- 



Now the MSSM Yukawa couphng constants are expressed as : 



2c. 



1 - 2c, 



k 



2{i - e-(i-2c.)fcR^} y 2{1 - e-(i-2c^)^«^} M5 

k 



1 - 2c. 



1 - 2c, 



d kl 



{ye)r 



2{l _ e-(i-2cfc)fci?,7r j Y 2{1 - e-{i-2c;)fcR'r j j^^ 

k 



2c 



2c 



2{1 - e-(l~2c„)fci?7r j y 2{1 _ e-(l-2c„)fcfl7r| j^^^ 

and the neutrino mass matrix nii, is given by : 



{yu)ij ; 

{yd)ki , 

{ye)mn 



1 -2c„ 



l-2c„ 



2{i _ e-(i-2cp)fci?7r} Y 2{1 - e"(i-2c9)fc«'r} 



v)pq 



(19) 



The geometrical factor ^^(1 — 2c) / (2{1 — e~(i-2c)fcR'r}) generates hierarchical couplings with- 
out fundamental hierarchy [2]; for c < 1/2, it is approximated by ^^1/2 — c and is 0(1), whereas 
for c > 1/2, it is approximated by \/ c — 1/2 e-(c"i/2)fc-Rvr jg exponentially suppressed. Note 
that this factor cannot be larger than 0(1). We assume that the components of 5D coupling 
matrices y„, y^, are all 0(1) and that the hierarchy of the Yukawa coupling constants stems 
solely from the following terms: 



2c,; 



29 



2{1 - e-(l-2ci)fcij7r j y 2{1 - e~(l~2Cj)fcR7r j " 

This is how the bulk matter RS model explains the Yukawa coupling hierarchy. 

We further assume that the components of Yy are 0(1). Note that Yy arises by integrating 
out singlet Majorana neutrinos. If the components of 5D neutrino Dirac coupling are 0(1), 
it is possible to take the value of Mgeesaw) such that {Yy)i^ ~ 0(1) holds, regardless of the 
5D disposition of singlet neutrino fields and the flavor structure of the Majorana mass term. 
Hence this is a natural assumption in the bulk matter RS model, in which all 5D couplings are 
considered 0(1). With this assumption, the hierarchy of the light neutrino mass matrix (19) 
arises only from the geometrical factors of SU(2) doublet lepton fields. 

We hereafter use the following notations: 



2c 



OLi 



qi 



2{1 _ f.-{l~2Cqi)kRfK^^ 



with i = 1, 2, 3 



(20) 
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for SU(2) doublet quark superfields with flavor index i, and 7^, Si, ej for SU(2) singlet up- 
type quark, singlet down-type quark, doublet lepton and singlet charged lepton, respectively. 
Then the hierarchical structures of the up-type quark Yukawa matrix Yu, the down-type Yukawa 
matrix Yd and the charged lepton Yukawa matrix Ye (in the basis of diagonal 5D bulk masses) 
are expressed as: 

~ Attj , (Xd)ij ~ 7j«j , (Xe)ij ~ ei6j , (21) 
and that of the neutrino mass matrix rrii, is expressed as: 

i^^h - ^r^3 M '\-^R^ ■ (22) 



2.2 Determination of the Geometrical Factors 

The magnitudes of the geometrical factors, a^, 7^, e^, can be almost determined by the 
data on SM fermion masses, CKM matrix and neutrino oscillations. The sole exception is 
the absolute scale of 6iS, of which we only know the relative scales between different flavors. 
In this subsection, we will estimate these factors. The values that correspond to the model 
must be given at the KK scale, Mse"^^'^, where the 5D theory is connected to the 4D effective 
theory. However, as is seen from [13], RG evolutions change the Yukawa coupling constants by 
at most 2 and the CKM matrix components by at most 1.2 through evolving from ~ 10^^ GeV 
to the electroweak scale. Also the neutrino mass matrix is affected only by 0(1) through RG 
evolutions [11]. Therefore we may estimate the magnitudes of at, Pi,'ji,6i,ei from the data at 
low energies. 

We first derive the model's predictions on the eigenvalues of the Yukawa coupling matrices 
and the components of CKM matrix. Let us diagonalize the Yukawa matrices: 

V^Y^Ul = diag, 
VdYdUl = diag, 
VeYMl = diag. 

For successful diagonalization of the hierarchical Yukawa matrices, the unitary matrices, 
Uu, Ud, Vu, Vd, Ue, Ve, need to have the following structure: 

/I 0\ 

Uu ^ Ud ^ «i/«2 1 , 14 ~ (a ^ /3) , ~ (a ^ 7) , 
\ ai/tts a2/tt3 1 / 

f/e ~ (a ^ 5) , V; ~ (« e) (23) 
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which leads to 



VuYuUl ~ diag ( /32a2, /^s^s ) , 

VdYdUl ~ diag ( 7iai, 72^2, Ts^s ) , 

VeYMl ~ diag ( ei5i, €2^2, £3^3 ) • (24) 
The hierarchical structure of CKM matrix Uckm is given by 

(1 ai/a2 oii/a-i \ 
axla-i 1 a2/«3 • (25) 
ai/as a2/«3 1 / 

We next list the experimental data on CKM matrix and the neutrino mass matrix. The 
absolute values of the CKM matrix components, \Uckm\-, at the electroweak scale has been 
measured to be [T5] 



\Vckm\Mw\\ 

0.97419 ±0.00022 0.2257 ± 0.0010 0.00359 ± 0.00016 

0.2256 ±0.0010 0.97334 ± 0.00023 0.0415 ± 0.0010 - 0.0011 

0.00874 ± 0.00026 - 0.00037 0.0407 ± 0.0010 0.999133 ± 0.000044 - 0.000043 

We approximate this matrix by the following formula: 

/I A A3 \ 

\Vckm\ ^ a 1 am with A = 0.22 . (26) 




To discuss the neutrino mass matrix, we adopt the tri-bi-maximal mixing matrix: 

/ ./i ./I n \ 

Umns = 



3 V3 



6 V 3 V 2 



V vi vi vi / 



and the following data on neutrino mass squared differences 

Aml-^ = 7.59 ± 0.20 x 10"^ eV^ \Aml^\ = 2.43 ± 0.13 x 10^^ gy^. 
We assume that the mass of the lightest neutrino is negligible. Then the neutrino mass matrix. 
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Umns diag ( m^i, m^2, fn^s ) ^Ijj^g , is given by 



/ 0.29 0.29 0.29 

Umns diag ( m,i, m,2, ) U^j^s = 0-29 2.8 -2.2 | x lO'^^ GeV 

\ 0.29 -2.2 2.8 

for normal hierarchy , 



4.9 0.026 0.026 

Umns diag ( m,i, m,2, m,-, ) Ulj,s = I 0.026 2.5 2.5 1 x 10"" GeV 



(27) 



0.026 2.5 2.5 



for inverted hierarchy . 



(2^ 



We now compare the predictions of the model with the data and estimate the magnitudes 
of Q!j, (3i, 'ji, 6i, Ei . For Yukawa eigenvalues, we simply have 

I3iai ~ m„/wsin/3, /32«2 ~ mc/t'sin/3 , P^a^ ~ m^/f sin/3 , (29) 
7iai ~ m^/f cos/3, 72^2 ~ ms/vcosf3, 73^3 ~ mb/vcos(3, (30) 
ei(5i ~ me/t'cos/3, 62^2 ~ m^/t>cos/3 , 63^3 ~ m^/t>cos/3, (31) 

where the mass values can be approximated by their pole values. Since the top Yukawa coupling 
is ~ 1, we have 03/33 ~ 1, which leads to 

03 ~ 1 , /33 ~ 1 . (32) 

Comparing (25) with (26), we find that putting 

ai ~ , ^2 ~ (33) 

works. We then have 

/3i ~ niu/v sin (3 , (32 ~ A^^ mc/f sin /3 , (34) 
7i ~ nid/ V cos (3 , 72 ~ A~^ m^/w cos/3 , 73 ~ m}y/vcos(3. (35) 

Next compare the matrix (22) with the neutrino mass matrix. For normal hierarchy case, it 
is possible to reproduce the hierachical structure of the neutrino mass matrix by assuming the 
relation: 

3(5i ~ ~ ^3 ; (36) 

and the ratio up to 3 among the components of 5D coupling Yy. In contrast, for inverted 
hierarchy case, ~ 200 ratio is required among the 5D coupling components no matter how 
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we choose SiS, which makes it difficult to naturally explain the hierarchy of the neutrino mass 
matrix. The situation gets worse if we consider non-negiligible mass of the hghtest neutrino. 
In conclusion, the bulk matter RS model favors the normal hierarchy of neutrino masses and 
imphes the relation (34) for 6iS. We estimate e^s assuming the relation (36); we obtain 

ei ~ 3 ^3'"'^ nie/v cos (3 , 62 ~ 5^^ m^/vcos(3 , €3 ~ 5^^ rrir/v cos (3 . (37) 

The magnitude of ^3 is a free parameter because we do not specify the scale of Mgeesaw 



3 Flavor- Violating Soft SUSY Breaking Terms 

In this model, flavor-conserving soft SUSY breaking terms arise from RG contributions of 
gaugino masses below the KK scale (gaugino mediation) and gauge interactions with messenger 
superfields (gauge mediation). On the other hand, flavor- violating terms arise from contact in- 
teractions with SUSY breaking sector on the IR brane (gravity mediation) and RG contributions 
of the Yukawa couplings. Of particular importance are the gravity mediation contributions, 
which have a flavor structure unique to the bulk matter RS model. In this section, we sepa- 
rately estimate the gravity mediation contributions and the Yukawa coupling contributions to 
flavor-violaing soft SUSY breaking terms. 

Remember that there are two scales of soft SUSY breaking terms, namely gravity mediation 
scale and gauge mediation scale. Assuming that the messenger masses and couplings are around 
the same orders, we define the following two mass parameters: 

1 Amp.^s I ^ Fj^ > 



M = --mess I ^ - A - I /oq^ 



where M^ess represents the typical scale of the SUSY conserving messenger masses Mmessii and 
Xmess the typical scale of the messenger couplings to SUSY breaking sector Xmessi- Note that 
Yukawa RG contributions to flavor- violating terms depend on both Mg^av and Mgaugei whereas 
gravity mediation contributions do only on Mg^av 



3.1 Flavor- Violating Gravity Mediation Contributions 

Let us estimate the magnitudes of gravity mediation contributions in the bulk matter RS 
model. 
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From (14), we obtain the following formulae for gravity-mediation-origined soft SUSY break- 
ing matter mass terms at the scale M^e'^^'^: 
For SU(2) doublet squarks, we have 

k 

i^Dij = i-dQ2ij + Jf^ ^grav ■ (40) 

By substituting {U, (3), {D, 7), (L, 6), {E, e) into {Q, a) in the above formula, we obtain similar 
expressions for SU(2) singlet up-type squarks, down-type squarks, SU(2) doublet sleptons and 
singlet charged sleptons. 

Assuming that the 5D couplings d^2ij, d^uj are 0(1), we obtain the following estimates on 
the magnitudes at the scale M^e~^^'^: 

{rriDij ~ aiaj M^,^^ . (41) 

and similar formulae with {U,P), {D,'~f), {L,6), {E,e) replacing (Q, a) in the above formula. 

Next we estimate the magnitudes of the A-terms that are induced by gravity mediation. 
The terms (15) directly contribute to the A-terms. Furthermore, since the Higgs superfields 
can couple to SUSY breaking sector in the following way: 

"^^^jrpnh^^nHu + ddAj^—^HlHd + h.c. 



the A-terms also arise from the Higgs F-terms via (12). Hence gravity-mediation-origined 
A-terms at the scale M^e~^^'^ are given by 

k k 
ik/5 ikzs 
k 

duA i^u)ij -^grav ~\~ iflu)ij f^i^^j , r -^grav i (42) 

k 

Adij = -ddA (Xd)ij Mgrav + {ad)ij littj — Mgrav , (43) 

M5 

k 

Aeij = -ddA {Ye)ij Mgrav + {ae)ij ei6j — Mgrav ■ (44) 

ik/5 

Assuming that the components of 5D couplings (i*A, {cL*)ij are 0(1), we obtain the following 
estimates on the magnitudes at the scale M^e~^^'^: 

Auij ~ {Yu)ij Mgrav + Pi<yj Mgrav , (45) 
Adij ~ {Yd)ij Mgrav + liOij Mgrav , (46) 
Aeij ~ {Ye)ij Mgrav + U^j Mgrav ■ (47) 
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3.2 RG Contributions 



Let us estimate the magnitudes of the flavor-violating soft SUSY breaking terms that arise 
from the RG equations involving the Yukawa couplings. In doing so, we take the specific flavor 
basis where or and Y^. are diagonal. 

We first study how y„, Y^., l^e-diagonal bases change through RG evolutions. Define the 
following unitary matrices f/*: 

UuYuUqu = {diag.) , 
UDYdUgd = (diag.) , 
UeYMl = (diag.) . 

Note that f/*s depend on the renormalization scale because the Yukawa matrices receive RG 
corrections. We will calculate how U^s vary through RG evolutions. We have the RG equation 
below: 

f-i^iUuY^Ug^) = {^A.u^)uI{UuYMqu) + Uu{fi-^Y^)UQ^ + {UuY^UQ.,)UlJfi^UQ^) 
= {f^-^Uu)UliUuYMQu) 

+ { "^uYjYd + SY^Yu + MYlYuK + tT[YlYD]Y^ 

+ iUuY^UQ^)UlJf^^UQu), (48) 

where Y^ is neutrino Dirac coupling which appears if there exist singlet neutrinos lighter than 
the KK scale, from which the RG equations are calculated. We hereafter adopt GUT normal- 
ization for gi. From (48), we see that, to keep UuYJJqu diagonal through RG evolutions, it is 
sufficient to have 

^i^Uu = 0, (49) 
d/i 

'"ci/I^'^" ^ ~ 1^ ^'^^ ~ components of Y^Yd) Uqu ■ (50) 
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In a similar manner, we obtain the following sufficient conditions for other U^s: 

^^Ud = 0, (51) 
d/i 

^J'-j-Uqd = -777-^ (off - diagonal components of FJF„) Ugd , (52) 

d/i lOTT^ 

fx^UE = 0, (53) 
d/i 

u-^Ul = 7 (off — diagonal components of yXYo) Ul . (54) 

d/i 167r^ 



Now that we know how Yu, Yd, Vg-diagonal bases change through RG evolutions, we estimate 
the RG contributions to the A-terms in these bases. Below is the list of the MSSM RG equations 
for the A-terms: 

WTr'fi^iUuA^Ugu) = 3UuA^Y^YMq^ + 3UuYXAuUqu 



d/i 



+ (f/f/A„f/Q J (diagonal part of U^^Y^YdUg^) + 2UuY^YjAdUQu 

+ 2( 3tr[FjA.] - ^glM^ff - ^alKn - y^?3K/f WuYuUq^) 

+ ( 3tr[F,ty^] _ 1^^2 _ 3^2 _ ^^2 )^UuA^Uq^) 

+ ii[YlYn\{UuA^UQ^) + ii[Y},An\{UuY.JjQ^) , (55) 



lQ'K^^^{UDAdUQd) = 3UDAdYjYdUQd + 3UDYdYjAdUQd 



d/j, 



+ (t/BArf?7Qrf) (diagonal part of U^^Y^YuUgd) + 2UDYdY^A^UQd 

7 16 
+ 2( 3tr[F;A,] + tr^tAJ - -gfM^ff " ^d^M^ff ' y^^s^vf WuYdJQd) 

7 16 

+ ( 3tr[riF,] + tr^tFe] - -gl - ^qI - -qI ){UDAdUQd) , (56) 



l&^l''^l^{UEA,UL) = 3UEA,Y}YJJL + 3UEYXAeUL 



d/j, 



+ 2( 3tT[YjAd] + tr[F;Ae] - ^^??Mf/f - 3glM^Jl^ ){UeYMl) 



+ ( 3tr[F]F,] + tT[YX] - - ^gl WeA^Ul) 

+ (t/^A^^i) (diagonal part of UIy},YdUl) + 2UeY,yIAdUl , (57) 

where Yd and Ae, respectively indicate neutrino Dirac coupling and its corresponding A-term. 
Note that the magnitudes of the components of the Yukawa couplings in each basis are 
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given by {Sij is the ordinary Kronecker's delta): 

{UuYJjQu)ij ~ Pitti 6ij , {UuYdUQu)ij ~ TjOj , 
{UDYuUQdlij ~ , {UDYdUQdlij ~ 7i«i f^ii ) 

{UEYeUL)ij ~ Cj^i 5ij , {UeYdUl)!] ~ O'^i 5 

where CjS indicate the geometrical factors for singlet neutrinos and satisfy 0^1- Note also 
that the A-terms receive RG corrections which are proportional to the corresponding Yukawa 
couplings and to the gaugino masses. We write these terms by M„, M^, Me, M/) respectively 
for Au, Ad, Ae, Ad- They depend on both Mgrav and Mgauge- 

With these ingredients, we estimate the RG contributions to those parts of A-terms which 
are not proportional to the corresponding Yukawa couplings, or equivalently the off-diagonal 
components of {UuAuUqu), (UDAdUQd), {UeAcUl)- In the right hand sides of (55-57), the 
second lines determine the magnitudes of the RG contributions. We thus obtain the following 
estimates (z 7^ j): 

A{UuAuUQuhj ~ 2 A(a.)'(73)'«, ^ ^ j ^(In/^) , (58) 

A{UDAdUQd)^, ~ 2 ii{aif{P^faj ^ ^ j d(ln/i) M„ , (59) 

/^{UEA,UL)^J ~ 2 e,(5,)2(C3)'5, X -i^ /d(ln/i)MB. (60) 

Since M„, M^, depend on Mg^av and Mgauge, so do the magnitudes of the RG contributions 
above. 



Let us estimate the RG contributions to soft SUSY breaking matter mass terms. Below are 
those parts of the MSSM RG equations that give rise to flavor-violating soft SUSY breaking 
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mass terms: 



WTi'fx—iUumljUl) 



+ (diagonal parts of Ul^YlYdUQu){Uljn%UQu) 
+ (diagonal parts of Ul^Y^YdUgu) 

+ 2Ul^YjmlY,UQ^ + 2iUlXY,UQ^)ml^ 

+ 2Ul^AiA^UQ^ + 2UI,4a,Uq^ , 



(61) 



(62) 



+ AUnY^mlYjuj, + AiUnY,Yiul,)ml^ 



(63) 



167r^fi—{UlmlUr_ 



+ 2UIy^Wj,Y,Ul + 2{UlYXUL)ml^ 

+ 2UIaIA,Ul 

+ (diagonal parts of UlY^YDUL){UlmlUL) 

+ (f/^m^f/i) (diagonal parts of UIy^Y^Ul) 

+ 2UlY^m%YDUL + 2(?7lF^F^f/L)m^^ + 2UlA^j,AnUL , (64) 



167rV^(t/£m|4) D 2f/^FeV';m|4 + 2C/^m|yeV;^[/; 



+ 4f/^yemiF;?7l; + AiUEY,Y^Uj,)ml^ 
+ AUeA^API,. 



(65) 



We first focus on the differences among the diagonal components of different flavors. From 
(61), the difference between the components {UQ^mQUQu)ii and {Uq^nQUQu)]] that arises 
through RG evolutions is given by (i > j) : 



A{ {Ul^u^lUQuh - {Ul^u^lUQ^h } ~ 2(a,)2(A)2 -L-^ J d(ln/i) ( 

+ 2(a,)2(73)2 -i-^ y d(ln/i) ( + m^, + mj,^ + M} 



(66) 
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where we neglected the terms proportional to {ajY because they are smaller than those propor- 
tional to (aj)^. Similarly, the difference between {U'lm^ULju and {,Ulm\UL)jj is given, from 
(64), by 



A{ {UlmlUL)u - {UWl^l),, } ~ 2 (5,)2(eO^ ^ | d(ln^) ( + m| + + M, 



+ 2 (5.)'(C3)' T7^ / d(ln/i) ( mi + + m^^ + 



167r^ 



(67) 



On the other hand, the differences among the diagonal components of SU(2) singlet soft mass 
terms follow different formulae. From (62, 63, 65), we have 

A{ {UjjmlUuh - (UlmlUu),, } ~ 4 (A)^(aO' J d(ln/x) ( + m| + m^^ + M,^ ) , 

(68) 

A{ (Uj^mlUnh ' (UlmlUnh } ~ 4 (7.)'(«i)' J d(ln/i) ( m^, + m| + m^,^ + Mj ) , 

(69) 

A{ (4m|f/^),, - iU],mlUE),, } ~ 4 (6,)2(5,)2 -i^ | d(ln/i) ( m| + + m^^ + ) . 

(70) 



We next study the off-diagonal components. In (61), terms 2UQji^Jm'j)YiiUQu, 2{UQ^YjYdUQu 
2UQ^A\AdUQu generate off-diagonal components, whose magnitudes are given by {i ^ j) 

A(f/(^,mJt/Qj., ~ 2 a,(73)'a, J d(ln/x) ( ml + + ) . (71) 
Similarly, we have 

AiUlmlULh ~ 2 5,(C3)'5, 1 d(ln/i) ( + m^^ + ) . (72) 

On the other hand, RG contributions to the off-diagonal components of {UjjmfjUu) arise from 
those of {Ui^^mlUQj and (UuA^Uqu) via terms 4?7c/r„m|yJf/J, AUuA^AlUl in (62). From 
(58), (71), (62), we obtain the following estimate on the magnitudes of the off-diagonal com- 
ponents {i 7^ j): 

A{UlmlUu)^, ~ 8 A(«.)'(73)'(«,)'/3. (j^) jdilnfi) J d{ln fi') { ml + mj,^ + M', ) 



]m 



+ 16 A(«.)'(73)'(«,)'/3j- (y^) y"d(ln/i) (^Muj d{\n ft')Md 



(73) 



1^ 



In the same way, we obtain the following estimates on the RG contributions to the off-diagonal 
components of (f/j^ml, (Ul^m'^UE)- 

AiUlmlUoh ~ 8 7.(a.)'(/33)'K)S (jL^) /d(ln^) J d{\n f,') { ml + m],^ + ) 
+ 16 7.(«.)'(/33)'K)S (y^) y"d(ln/i) (^M, I di\n fi')Mu^ , (74) 

A(f/l;m|^^),, ~ 8 e.(5.)2(C3)'(5,)%- (j^) y"d(ln/.) J d{\n ^^') { m% + mj,^ + M^) 
+ 16 e,(5,)2(C3)^(5,)% y"d(ln/i) (^M, j d{lnf,')Mo^ . (75) 



Finally, we briefly discuss whether this model gives a realistic mass spectrum consistent 
with the bounds on flavor-violating processes. 

For cases without messenger flelds, i.e. when gaugino mediation is the only source of soft 
SUSY breaking masses, the paper [5] showed that there exist mass spectra below TeV scale 
that satisfy all experimental bounds. However, ~ 0.1 suppression on the term {Ag)2i relative 
to its natural scale (~ G2^iMgrav) is required to evade the bound on /i — > 67 branching ratio. 
Other soft SUSY breaking terms are less constrained. 

If there are 1 to several messenger pairs, the resultant mass spectra are more likely to evade 
the experimental bounds because gauge mediation contributes solely to flavor-universal soft 
SUSY breaking terms. 



4 Signatures of the Model 

In the previous section, we saw that the bulk matter RS model combined with 5D MSSM 
predicts a unique flavor structure of gravity mediation contributions to flavor-violating soft 
terms. We here discuss the ways to observe this structure through future collider experiments. 

Focus on the flavor compositions of SUSY matter particle mass eigenstates. Due to flavor- 
violating soft mass terms {ml)ij and flavor-violating A-terms (A*)jj, SUSY particles of different 
flavors mix in one mass eigenstate, whose flavor composition reflects the relative size of the 
flavor- violating terms. Since sparticles of different flavors decay into different SM particles (plus 
the lightest or the next-to-lightest SUSY particle), one can measure the flavor composition by 
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detecting the decay products of that mass eigenstate, counting the event numbers of different 
decay modes and calculating their ratios. These ratios are connected to the structure of flavor- 
violating soft SUSY breaking terms and make it possible to experimentally test the predictions 
of the bulk matter RS model. 

Below we formulate the relation between flavor-violating terms and sparticle flavor mixings. 
In the first subsection, we interpret the predictions of the bulk matter RS model in terms of 
the flavor mixing ratios of sparticle mass eigenstates. In the next subsection, we look into the 
predictions of models other than the bulk matter RS model and discuss whether or not it is 
possible to distinguish different models. 

Consider the situation where sparticle "a" with soft SUSY breaking mass mixes with 
sparticle "b" with soft mass ml through mixing term Am^. The mass matrix in the basis of 
(a, b) is given by 

/ ml Am^ \ 
\^ Am^ ml J 

The mass eigenstates are derived by diagonalizing the matrix above. If \ml — ml\ » 2|Am^| 
holds, the mixing ratios of "a" and "b" in the two mass eigenstates are approximately given by 

\ml — ml\ : |Am^| , |Am^| : \ml — ml\ . 



4.1 Predictions of the Bulk Matter RS Model 

The bulk matter RS model predicts a nontrivial structure of flavor- violating soft SUSY 
breaking terms, given by (41, 45-47, 66-75). This structure can be translated into the flavor 
composition of each SUSY particle mass eigenstate. One subtlety is that the flavor-violating 
terms contain two different SUSY breaking mass scales, namely, the IR-scale-suppressed gravity 
mediation scale, Mgrav, and the gauge mediation scale, Mgauge] flavor-violating gravity media- 
tion contributions depend solely on Mg^av, whereas RG contributions are proportional to the 
net soft SUSY breaking mass scale that depends both on Mg^av and M gauge- The relative size 
of these scales affects the predictions on the flavor compositions. We consider three cases with 
Mgrav > M gauge, Mgrav < M gauge and Mgrav « M gauge, whose precise definitions wiU be given 
each time. These cases lead to different predictions. 
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4.1.1 Case with Mgrav > M gauge 

In this case, flavor-universal soft SUSY breaking masses, m^, and gaugino masses, M^,*, are of 
the same magnitude as the gravity mediation scale, Mg^av The differences among the diagonal 
components of different flavors come from the gravity mediation contributions (41) and the RG 
contributions (66-70). Since we now have ml ~ Mg^.^^, M*^, ~ Mgrav, terms (41) surpass terms 
(66-70). Hence we may make the following approximations for i > j in any flavor basis: 

(mj)ii - {ml)jj ~ M^^av (76) 

and similar formulae with {U,/3), {D,'-/), {L,6), {E,e) relacing {Q, a) in the above formula. 
In a similar manner, in any flavor basis, the A-terms are approximated by 

{Au)ij D Attj Mgrav , {Ad)ij 3 ^iOj Mgrav , (^e)ij 3 ei6j Mgrav , (77) 

and the off-diagonal components of soft SUSY breaking mass terms are by {i ^ j) 



/I] Lxjuj ^"---grav 



and similar formulae with {U,I3), {D,'~f), {L,6), {E,e) relacing (Q, a) in the above formula. 

Sparticle Qi mixes with sparticle Qj {j ^ i) through the term {rn^)ij and with Uk or 
Dk {k 7^ i) through the A-terms and the VEVs of the Higgs bosons. In this way, there appears 
a mass eigenstate that consists mainly of Qi and partly of Qj and Uk or Dk, which we hereafter 
call "almost Qi mass eigenstate". From (76, 78), the mixing ratio of Qj in "almost Qi" mass 
eigenstate is given by 



I("^q)uI (^if^jM^rav _ aj ^^g^ 



for i > j, and by 



{ml)ii- {ml)jj\ (tti)^^gm,; " 

\iml)ij\ aiUjMlrav 



~ — (80) 
a,- 



I {ml)ii - {ml)jj I {ajfM'^rav 

for i < j. On the other hand, the mixing ratio of Uj in the up-sector of "almost Qi mass 
eigenstate" is given by (z 7^ j) 

|(^«)ji| f^j^^iM^grav 



T7^ , (81) 



where we used the fact that the difference between the flavor-universal masses of SU(2) doublet 
and singlet squarks is of the same magnitude as the soft SUSY breaking mass scale itself, 
denoted by Mgusy 



21 



The mixing ratios in other mass eigenstates follow similar formulae. There is a subtlety 
about the ratio of Lj in "almost Lj mass eigenstate" because we have 3(5i ~ ^2 ~ ^3 and the 
approximation used to derive (79-81) is no longer valid. Actually, the mixing ratio of Lj in 
"almost Li mass eigenstate" is 0(1) for any z, j. 



4.1.2 Case with Mgrav « M gauge 

In this sub-subsection, we concentrate on the case where the ratio Mg^av/M gauge is so small 
that the RG contributions to flavor- violating soft SUSY breaking terms, (58-60, 66-75), are of 
the same magnitude as or larger than the gravity mediation contributions, (41, 45-47). 

In these cases, the mixing ratio of Qj in "almost Qi mass eigenstate" is given, from (66, 
71), by {i>j) 

aiilzfaj «i (73)' .^21 

IK)..-Kkl ^ {a^m^Y + {a,ni,f ^ a. (A)^ + (73)^ ' ^ ^ 

in the flavor basis where Yu is diagonalized. Here we used the fact that the integrand of the 
right hand side of (66) and that of (71) are of the same magnitude. On the other hand, from 
(58), the mixing ratio of Uj in the up-sector of "almost Qi mass eigenstate" is given by {i 7^ j) 

1^2 _^2| ~ 2/3,(a,) (73) a, — (83) 

\ Q '"'Ul ^''J-gauge 

in y^-diagonal basis. Here we approximated the difference of flavor-conserving masses of SU(2) 
doublet squarks and singlet up- type squarks by M gauge- The mixing ratio of Dj in the down- 
sector of "almost Qi mass eigenstate" in y^-diagonal basis takes a similar expression. The same 
discussion applies to the mixings in "almost Lj mass eigenstate" . 

The mixing ratios in "almost Ui mass eigenstate" follow different formulae. From (68, 73), 
the ratio of Uj is given by {i > j) 

\{ml),^\ 24 A(a,)2(73)2(a,)2/3^- , .,(3^ 



6 (73)^(«.)^ ^ (84) 



\{mlh-iml),,\ 4(A)2(«,)2 ^'^^ ^ A 

in y„-diagonal basis. On the other hand, from (58), the ratio of the up-sector of Qj in "almost 
Ui mass eigenstate" is given by {i 7^ j) 

1^2 _^2| ~ 2A(«^) (73) aj — (85) 

\"''Q '"'Ul ^"gauge 

in y^-diagonal basis. The same discussion applies to the mixings in "almost Di mass eigenstate" 
and "almost Ei mass eigenstate". 
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4.1.3 Case with Mgrav < M gauge but not with Mg^av « M gauge 

Consider the case where Mg^av is slightly smaller than M gauge- Then gravity mediation 
contributions surpass RG contributions for some of the flavor-violating soft SUSY breaking 
terms, and the opposite holds for the other terms. In these cases, the mixing ratios of sparticle 
mass eigenstates generally depend on the unknown ratio Mgrav /Mgauge and the model loses its 
predictive power. 

However, certain mixing ratios are more likely to reflect the gravity mediation contribu- 
tions. For example, if Mgrav ^ ^zMgauge, as to terms {rn\)ii — {rny)jj and {mE)fj, the 
gravity mediation contributions described by (41) are larger than the RG contributions, (70, 
75). Then the mixing ratio of Ej in "almost Ei mass eigenstate" is the same as in the case with 
Mgrav ^ Mgauge- Focusiug ou such mixing ratios, it is still possible to observe the signatures of 
the model. 



4.2 Comparison with Other Models 

To test the predictions of the bulk matter RS model, we must check whether they con- 
tain signatures distinguishable from other models. As an example, we investigate two types of 
models; one is "minimal flavor violation" , in which RG contributions of the Yukawa couplings 
are the only source of flavor-violating soft SUSY breaking terms. The other is "4D gravity 
mediation" , in which gravity mediation contributes uniformly to all flavor- violating terms. We 
will compare the predictions of these models with the bulk matter RS model and discuss the 
ways to distinguish them. 



4.2.1 Minimal Flavor Violation 

The minimal flavor violation (MFV) scenario leads to the same result as in section 4.1.2, i.e. 
the bulk matter RS model with Mgrav << Mgauge- This is because the argument in section 4.1.2 
holds irrespective of gravity mediation contributions. We thus conclude that it is impossible 
to experimentally distinguish the bulk matter RS model from the MFV scenario when we have 

Mgrav « Mgauge, aS lu 4.1.2. 

In contrast, if Mgrav ^ Mgauge, the MFV scenario and the bulk matter RS model have 
distinctively different predictions on the mixing ratios in "almost Ui, Di, Ei mass eigenstates" 
with 2 = 1,2. This is seen by comparing (79-81) [Q relaced by [/, -D, E) with (84, 85); the flavor 
mixings in these mass eigenstates are suppressed at least by («2)^ or (^2)^ in "minimal flavor 
violation" compared to the bulk matter RS model. Therefore it is possible to discriminate the 
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two models by observing the flavor compositions of "almost 1st or 2nd generation SU(2) singlet 
sparticle mass eigenstates" . 

4.2.2 4D Gravity Mediation 

We here discuss the case where 4D theory description is valid even at the Planck scale, or all 
matter superfields are confined on the same 4D brane. Then the gravity mediation contribu- 
tions are of the same magnitude irrespective of flavors. Of particular interest is the situation 
where the gravity mediation contributions surpass the flavor-violating RG contributions, which 
is the case when Mgrav is only slightly smaller than M gauge- In this situation, the differences 
between diagonal components of soft SUSY breaking masses {mDa — {rriDjj, and off-diagonal 
components (m^)jj, in any flavor basis are of the same magnitude. Then the mixing ratios in 
sparticle mass eigenstates are all 0(1). It is easy to distinguish this model from the bulk matter 
RS model, where the mixing ratios of recessive flavors are suppressed by the geometrical factors. 

5 Experimental Studies 

In the previous section, we saw that the bulk matter RS model has a unique prediction on 
the flavor compositions of sparticle mass eigenstates that may be distinguishable from other 
models. In this section, we study how to measure the predicted mixing ratios through collider 
experiments. We focus on the case where Mgrav — Mgauge holds or Mgrav is slightly smaller 
than Mgauge, and put emphasis on distinguishing the bulk matter RS model from the MFV 
scenario. 

The basic strategy is to create a specific mass eigenstate(s) of SUSY matter particles, detect 
its decay products and count the numbers of events of different decay modes. The branching 
fractions of different modes reflect the flavor composition of that mass eigenstate. There are, 
however, three challenges for this study. 

First, we have to detect small flavor components of sparticle mass eigenstates, which means 
that we have to observe rare decay events in collider experiments. For this purpose, the 
probability of misidentifying the decay products of the dominant mode as those of a rare 
mode must be negligibly small. For example, the stau components of "almost smuon mass 
eigenstates" are detectable because SM tau from the stau components, when we focus on its 
hadronic decay, leaves a signal different from muon events. However, it is impossible to observe 
the smuon components of "almost stau mass eigenstates" because SM tau from the dominant 
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stau components may decay into SM muon, which mimics the smuon component signaL 

Second, we have to extract the decay products of a specific mass eigenstate in order to 
compare the data with the predictions of the bulk matter RS modeL It is thus required to 
produce only specific mass eigenstates at a collider. This is achieved by lepton colliders, such 
as the ILC [16] and the CLIC [17], where the center-of-mass energy of a process is fixed. For 
example, the flavor-mixing ratios in "almost SU(2) doublet smuon mass eigenstate" and in 
"almost SU(2) singlet smuon mass eigenstate" are predicted to be different. To confirm this 
prediction, we must produce one of the two eigenstates selectively. If the latter is lighter than 
the former, we take the center-of-mass energy between their thresholds so that only the latter is 
created on-shell. We then measure the mixing ratios of the latter eigenstate through its decay 
products. In conclusion, lepton colliders are essential when studying the flavor compositions of 
sparticle mass eigenstates. 

Finally, we have to focus on "almost SU(2) singlet mass eigenstates" in order to discriminate 
the bulk matter RS model from the MFV scenario. This is understood by comparing the 
predictions of the bulk matter RS model, (79-81), with those of the MFV scenario, (82-85). 
Remember that we have 

73 ~ tan/3 — , /^s ~ 1 

V 

and we do not know the magnitude of ^3. Hence it can be the case that the mixing ratios of 
Qj in "almost Qi mass eigenstate" , and those of Lj in "almost Li mass eigenstate" are of the 
same magnitudes for the bulk matter RS model and the MFV scenario. In contrast, the mixing 
ratios of Uj,Dj,Ej in "almost Ui,Di,Ei mass eigenstates" are of the different magnitudes for 
the two models because the mixing ratios in the MFV scenario, (82-85), are suppressed by the 
factors (tti)^, («2)^, (<^i)^ or (^2)^ compared to those in the bulk matter RS model, (79-81). We 
further notice that "almost 3rd generation sparticle mass eigenstates" are not suitable for our 
study because the 3rd generation sparticles have significant left-right mixing terms due to their 
large Yukawa couplings. We conclude that observing the rare decays of "almost SU(2) singlet 
1st and 2nd generation mass eigenstates" is the only way to distinguish the bulk matter RS 
model and the MFV scenario. 

Taking these points into account, we will discuss three types of experiments that are feasible 
at future lepton colliders. The first type of experiment deals with the rare decay of "almost 
SU(2) singlet smuon mass eigenstate" into SM tau, which reflects the mixing of singlet smuon 
with stau. Another type of experiment deals with the rare decay of "almost SU(2) singlet smuon 
mass eigenstate" into SM electron or that of "almost SU(2) singlet selectron" into SM muon, 
which reflects the mixing of singlet smuon and selectron. The other type of experiment deals 
with the rare decay of "almost SU(2) singlet scharm mass eigenstate" into SM top, which reflects 
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the mixing of singlet scharm with stop. For a concrete discussion, we make assumptions on the 
SUSY particle mass spectrum in section 5.1. We then look into the three types of experiments 
in section 5.2., 5.3. and 5.4. 

"Almost the lighter stau / stop mass eigenstate" are schematically denoted by fi / ti, and 
"almost singlet selectron / smuon / scharm mass eigenstate" are by / jiR / jiR. 

It is impossible to do these experiments at hadron colliders. This is fundamentally because 
we need to create "almost SU(2) singlet mass eigenstates" exclusively, without creating their 
"almost SU(2) doublet" counterparts, in order to discriminate the bulk matter RS model from 
the MFV scenario. Hadron colliders would necessarily create both eigenstates, and the de- 
cay products of the latter would contaminate the signals that allow us to distinguish the two 
scenarios. It is true that "almost SU(2) singlet eigenstates" are normally lighter than their 
"almost SU(2) doublet" counterparts, and hence the production cross sections of the latter are 
lower even at hadron colliders. However, since the two scenarios only predict the orders of 
magnitudes of the branching ratios of rare events, even small contamination from the latter 
would make it difficult to test the predictions. 

5.1 Assumptions on the Mass Spectrum 

In this subsection, we make plausible assumptions on the SUSY particle mass spectrum that 
are consistent with the bulk matter RS model combined with 5D MSSM. 

We assume that squarks are heavier than sleptons and gluino is heavier than Wino and 
Bino because of their SU(3) charges. Also, Wino is assumed heavier than Bino due to its SU(2) 
charge. SU(2) doublet squarks are heavier than singlet squarks, and doublet sleptons are than 
singlet sleptons. Since gauge superfields are fiat in the bulk, i.e. they have no dependence, 
they obtain large soft SUSY breaking masses through contact terms on the IR brane. Therefore 
gluino tends to be heavier than squarks. Wino and Bino are heavier than sleptons but lighter 
than squarks. 

The /i-term is assumed larger than Wino and Bino masses, as is normally the case. 

We do not specify the mass order among doublet and singlet squarks and sleptons because 
gravity mediation contributions may distort the mass spectrum. However, we expect that 
the masses of the 1st and 2nd generation SUSY particles are almost degenerate because their 
Yukawa couplings as well as their overlap with the IR brane are small. 

Gravitino is always the lightest SUSY particle (LSP) because its mass is given by ~ 
TeV xe-^^'^. The next-to-lightest SUSY particle (NLSP) is "almost SU(2) singlet selectron", 
"almost singlet smuon" or "almost the lighter stau" mass eigenstate. The lifetime of the NLSP 
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satisfies 

ctNLSP - 487r ^ ~ 487r^^- 

m^LSP [mNLSPr 

We assume that the hfetime is enough long that the NLSP reaches the inner detector before it 
decays. 

The order of the sparticle soft SUSY breaking masses is summarized below: 

Hu, Hd > g > Ql > Qr > Xi, xl (■=. W) > Xi (■=. B) > II > Ir > ^3/2 • 



5.2 Type I - Smuon Rare Decay with Stau-like NLSP 

Consider the case where "almost the lighter stau mass eigenstate" (fi) is the NLSP. Tune the 
center-of-mass energy of the lepton collider between the thresholds of "almost singlet selectron 
/ smuon mass eigenstates" (eij//i/j) and "almost SU(2) doublet selectron / smuon", "almost 
the heavier stau" mass eigenstates. Then cr, (Ir, fi are produced on-shell, while other sparticle 
mass eigenstates are not. 

The signal for cr or fiR pair production is a pair of long-lived charged massive particles, 
which are NLSP fis, plus two pairs of hard SM leptons. Note that, since the masses of cr and 
fiR are almost degenerate, we cannot detect SM leptons emitted when the heavier one decays 
into the lighter one. Normally, we have two SM muons or electrons plus two SM taus in these 
events (we call this "main mode"), e.g. 

ee — JXr iIr — )■ fi t fi fi t fi (86) 

for the sumon production. However, due to the small stau components in cr / fxR, we may also 
have one SM muon or electron plus three SM taus in these events (we call this "rare mode"), 
e.g. 

ee — )■ fiR fiR —7- r r fi /i r f 1 (87) 

for the smuon production. Requiring hadronic decay of SM taus and taking advantage of the 
tau vertexing, one can reduce the probability of misidentifying a main mode event as a rare 
mode event to a negligible level. 
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The branching ratio of the rare mode is proportional to the square of the mixing ratio. 
The stau component in fiR plays a dominant role because the stau component in is much 
more suppressed. From (79) with {Q,a) replaced by {E,e), and from (81) with {Q,U,a, /3,Vu) 
replaced by {L, E,6,e,Vd), the bulk matter RS model predicts that the branching ratio of the 
rare mode is given by 

Br{flR T T fi) ~ ( — ) + ( €26, ^'^ 



•-susy 



J \-^susy 

where the first term comes from the mixing with singlet stau and the second from the mixing 
with doublet stau. If en is lighter than fiR, the branching ratio is reduced by 1/2 compared to 
the opposite case because a half of fl^s decay into e/js. However, this does not affect the order 
estimate above. Since we have Mg^sy » the second term is negligible and the branching 
ratio becomes 

Br{fiR ^ T T fi) ~ 0.004 . (89) 

Note that the prediction above may change by 0(0. 1) — 0(10) because we only know the 
magnitudes of the higher-dimensional couplings for soft SUSY breaking terms. 

On the other hand, the MFV scenario predicts that the branching ratio of the rare mode is 
given by 



BrifiR ^ rrn) ~ jQ(C,y(S,y^). + 2(C3)'(52)'e253 



~ ((3^2)" X 0.1 , (90) 

where we used (84) and (85) with (Q, U, a, /3, 7, v^) replaced by (L, E, 5, e, C, Vd). Although we 
cannot determine the magnitude of C3^2) we expect it to be smaller than 0.1; from (72), we 
have the following flavor- mixing term for SU(2) doublet smuon and selectron: 

(mi)i2 ~ 5i(C3)'52 ML, ~ ^((3^2)' ML, . 

For example, with M^usy = 500 GeV, m? = 500 GeV, = = 750 GeV, = 1000 
GeV, tan/3 = 10 and vanshing A-terms, taking (-^62 = 0.1 would saturate the current bound 
on yU — )■ 67 branching ratio, Br{fi — 67) < 1.2 x 10^^^ [18]. Hence the branching ratio of the 
rare mode in the MFV scenario satisfies 

Br{jlR ^ rrn) < 10"^ , (91) 
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which may get smaller if the bound on Br{fi — t- 67) improves. We conclude that the branching 
ratio of the rare mode is distinctively smaller in the MFV scenario than in the bulk matter RS 
model. 



5.3 Type II - NLSP Selectron Rare Decay into Muon, or NLSP 
Smuon Rare Decay into Electron / Tau 

Consider the case where "almost singlet smuon" or "almost singlet selectron2 mass eigenstate 
(/iij or cr) is the NLSP and is long-lived. Tune the center-of-mass energy slightly above the 
threshold of flji / cr so that they are produced with a low (5 (Lorentz velocity). Slow long-lived 
sleptons may be trapped in the inner detector. According to the paper [19], taking (3 < 0.2 
is sufficient to trap 600 GeV or lighter sleptons in the inner detector. We analyze the decay 
products of these sleptons to study their flavor compositions. 



First study the case where cr is lighter than jiR and is the long-lived NLSP. cr mainly decays 
into a SM electron and a gravitino (main mode). However, due to its smuon component, it also 
decays into a SM muon and a gravitino (rare mode). Hence we expect to observe rare mode 
events where one of the sparticle pair produced by the collider decays into a SM muon and the 
other into a SM electron with large vertex separation due to the longevity of cr. 

The bulk matter RS model predicts that the branching ratio of the rare mode is given by 



rrie \ ^ (nip 
3— + 3 



'^/j / V -^susy 



0.0002 , (92) 



where we neglected the second terms of the right hand sides because we have 



3^ » 3- 



'^/i -^susy 

in realistic SUSY models. 

On the other hand, the MFV scenario predicts that the branching ratio of the rare mode is 
given by 

Br{eR ^ ^^ ~ |6(C3)'(5i)'-| + [2{C,f{5,fe,5,^\ 

L ^2 ) \^ susy ) 

~ (Cs^i)' X 0.03 , (93) 
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Again, the bound on /i — > 67 branching ratio gives a severe constraint on the value of Cs^i; cind 
the branching ratio satisfies 

Br{eii ^ /i^3/2) < 10-' (94) 

for reahstic mass spectra. 

cr also decays into a SM tau and a gravitino but the branching ratio is suppressed by the 
factor (mg/mT-)^ and is thus negligibly small. 

Next consider the case where /i/j is lighter than and is the long-lived NLSP. fiR mainly 
decays into a SM muon and a gravitino (main mode), but also into a SM electron and a grav- 
itino, or into a SM tau and a gravitino (rare modes). The branching ratio of the rare mode 
where the sparticle pair decay into a muon and an electron and two gravitinos is the same as 
(92, 93). The branching ratio of the rare mode where the sparticle pair decay into a tau and a 
muon and two gravitinos is the same as (89, 90). 

5.4 Type III - Scharm Rare Decay into SM Top 
5.4.1 Scharm is Lighter than Stop 

Consider the case where cr is lighter than ti. Tune the center-of-mass energy between the 
thresholds of cr and ii. Then cr, ur, sr, dR, whose masses are almost degenerate, are produced 
on-shell, while other squark mass eigenstates are not. cr / ur / sr / (Ir mainly decay into SM 
charm / up / strange / down and the lightest neutralino Xi? which is Bino-like (main mode). 
Xi promptly decays into several SM leptons and NLSP, e.g. we have 

ee ^ crCr ^ cxicxi ^ (c-jet) NLSP (c-jet) NLSP (SM leptons) (95) 

for scharm pair production. Due to the small stop components, they also decay into SM top 
and Xi with a tiny branching ratio (rare mode), e.g. we have 

ee crCr tXiCXi ^ (top decay products) NLSP (c-jet) NLSP (SM leptons) (96) 

for scharm pair production. The signal of the main mode is two hard jets, two long-lived charged 
massive particles and several SM leptons. On the other hand, the signal of the rare mode is, 
when SM top decays hadronically, four hard jets, two long-lived charged massive particles and 
several SM leptons. We see that the probability of misidentifying a main mode event as a rare 
mode event is negligibly small if we require the hadronic top decay in rare mode events. 
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Of the four eigenstates, cr dominantly contributes to rare mode events because the stop 
components in the other eigenstates are more suppressed than in cr. On can confirm this by 
requiring c-tagging for one of the jets in rare mode events. 

From (80) with {Q,a) replaced by {U, /3) and (81), the bulk matter RS model predicts that 
the branching ratio of the rare mode is given by 



Br(cR ^ t X?) ~ ( ^ ) + ( /32«3 



) V Msusy 

1 \ 2 /I 

1 nic \ I I nic rrit 



A2 nit J rrit Msusy, 

~ 0.02 , (97) 
where we neglected the second term because we have 

rrit < Msusy 

in realistic SUSY models. 

From (84) and (85), the MFV scenario predicts that the branching ratio of the rare mode 
is given by 

Bridn^tx'i) ~ {6(73)'(«2)'|} + {2(73)'(«2)'/32«3^} 

< 5 X 10-6 , (98) 

where we used the fact that 73 < 1. Comparing (96) with (95), we notice that the bulk matter 
RS model and the MFV scenario have distinctively different predictions. 



5.4.2 Stop is Lighter 

Consider the case where ti is lighter than c/j. cr is still lighter than "almost SU(2) dou- 
blet squark" mass eigenstates (including ti,6i). Tune the center-of-mass energy between the 
thresholds of cr and "almost doublet squark" eigenstates. Then cji,UR,Sfi,dR,bi,ii are pro- 
duced on-shell, while the other squark mass eigenstates are not. 

We want to extract the signals of rare mode events where one of the pair of cr, ur, sr or 
(Irs decays into a SM top and a neutralino. However these events are contaminated by the 
events where one of the pair of bis decays into a SM top and a chargino; the chargino decays 
into a NLSP and SM leptons, but one charged lepton is mis-detected. &-jet from the other bi is 
mis-6-tagged. There is also a contamination from the events where one of the pair of tis decays 
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into a SM bottom and a chargino, or into a SM charm and a neutralino due to the scharm 
component in ti. We take advantage of kinematical properties to reject these contaminations. 

Tune the center-of-mass energy close to the threshold of cr, mr, s/j, (i/j so that they are 
produced almost at rest. Suppose that one observed three hard jets from the hadronic decay of 
SM top (t), another hard jet (j) and several SM leptons plus two NLSPs. Further assume that 
the 3-momenta of t and j are reconstructed successfully. We want to know whether this event 
comes from the decay of C/j, Ur, Sr, dji or from ti, 61. In the former case, the 3-momenta satisfy 

\Pj\ + \/\ -Pj\^ + m?^ ^ rn-cji , (99) 
\l\pt? + ml + y/| -pi|2 + m2 - mg^ , (100) 

where pj and pt respectively denote the 3-momenta of j and t, and the mass of Bino-like 
neutralino. In the latter case, however, the above equations hold for specific situations where 
j and t go in special directions against the initial tis or feis because the ti, feis are boosted. To 
summarize, the rare mode signals of Cij, mr, sr^ dn can be extracted through the discriminants 
(99, 100). 

The branching ratios of the rare mode in the bulk matter RS model and the MFV scenario 
are the same as in (97, 98). 

5.5 Cross Sections 

We calculate the cross sections of the rare modes given above and discuss their accessibility 
at collider experiments. 

First focus on the Type I experiment, where one of the pair of "almost singlet smuon or 
selectron mass eigenstates" decays into two SM taus and a NLSP in case the NLSP is stau-like. 
The center-of-mass energy is tuned above the threshold of fiR, namely we take 

= 2m^^ + 100 GeV . 

In Figure 1, we plot the mass of ftR vs. the cross section of the rare mode at a e+e~ collider. 
We take the branching ratio of the rare mode as 



Br{fiR ^ r T Ti) = 




based on (88). We assume that fiR is slightly lighter than e/j so that a half of e/js decay into 
fiR and contribute to the rare mode. In calculating the cr pair production cross section, the 
Bino mass is assumed to be 1.5m^^. Also shown is the total cross section of flu, cr production 
events. 
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Figure 1: Mass of vs. the cross section of the rare mode where a fiR is produced and decays 
into two SM taus and a NLSP stau (the straight hue). The center-of-mass energy is taken as 
y/s = ^nifi^ + 100 GeV. The total cross section of fiR or cr production process is also shown 
(the dashed line). 

Next focus on the Type II, where one of the pair of "almost singlet selectron mass eigen- 
states" decays into a SM muon and a gravitino if it is the NLSP, or one of the pair of "almost 
singlet smuon mass eigenstates" decays into a SM electron / tau and a gravitino if it is the 
NLSP. We tune the center-of-mass energy slightly above the threshold of cr / jiR so that their 
velocities are low enough to trap them inside the inner detector. For simplicity, we take 



In Figure 2, we plot the mass of cr vs. the cross section of the rare mode if it is the NLSP, and 
the mass of jiR vs. the cross section of the rare mode if it is the NLSP. We take the branching 
ratios of the rare modes as 



based on (88, 92). Also shown are the total cross sections of cr and jiR production processes. 

Finally focus on the Type III, where one of the pair of cr, ur, sr, cIrs decays into a SM top 
and a neutralino. The center-of-mass energy is tuned above the threshold of cr. First we take 

y/s = 2me^ + 100 GeV 

so that the cross section is nearly maximized. Second we take 

y/s = 2mc^ + 10 GeV 



2"^efl/Afl + 20 GeV , 



Br{eR /i ^3/2) 
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Figure 2: Left: Mass of cr vs. the cross section of the rare mode where a NLSP cr is produced 
and decays into a SM muon and a gravitino (the straight hne). Right: Mass of flu vs. the 
cross section of a rare mode where a NLSP fxR is produced and decays into a SM electron and 
a gravitino (the straight hne below), and the other rare mode where a NLSP flu is produced 
and decays into a SM tau and a gravitino (the straight line above). The total cross sections of 
ept (left) and flji (right) production processes are also shown (the dashed lines). 

so that crs are produced almost at rest and the rare mode events are kinematically distinguish- 
able from ti, bi production events in case ti, bi are lighter. In Figure 3, we plot the mass of cr 
vs. the cross section of the rare mode for both cases. The branching ratio of the rare mode is 
taken as 

/I \ 2 
^ /- riN I 1 nir ^ 

Br{cR ^ txi) = 



rrit 

for both cases, based on (97). Also shown is the total cross section of cr, ur, sr, (Ir production 
processes. 
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Figure 3: Mass of cr vs. the cross section of the rare mode where a cr is produced and 
decays into a SM top and a neutralino (the straight lines). The center-of-mass energy is taken 
as ^/s = 2mcj^ + 100 GeV in the left and i/i = 2mc^ + 10 GeV in the right. The total cross 
section of cr, ur, sr or dR production processes is also shown (the dashed lines). 



Note that we can in principle reject all background events when observing the signals of 
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the rare modes. Therefore detecting several signals is sufficient to con&m the bulk matter RS 
model. From Figure 1 and 2, we find that one can study the stau component in /Tr at the ILC 
with the integrated luminosity of ~ lOOfb"^. However, studying the smuon component in 
or the selectron component in fipi requires ~ lOOOfb"^ integrated luminosity. From Figure 3, 
we see that ~ lOOfb"^ integrated luminosity is sufficient to study the stop component in c/j. 

6 Summary and Outlook 

We discussed observing signals of the bulk matter RS model, especially when the IR scale is 
far above TeV scale. We saw that this is possible in the case of the minimal supersymmteric 
extension of the bulk matter RS model where the warped spacetime solely explains the hierarchy 
of the Yukawa couplings, while SUSY solves the gauge hierarchy problem. There, gravity 
mediation contributions to soft SUSY breaking terms refiect the 5D disposition of superfields. 
Hence fiavor-violating soft SUSY breaking matter mass terms that arise from gravity mediation 
exhibit a fiavor structure unique to the bulk matter RS model. RG running of the Yukawa 
coupling constants also contributes to the fiavor-violating terms, but its contributions and the 
gravity mediation contributions are distinguishable if the mass scale of gauge mediation is not 
much larger than that of gravity mediation. Then the latter contributions can be extracted 
by investigating the 1st and 2nd generation SU(2) singlet SUSY particles, where the former 
are further suppressed by the small Yukawa coupling constants. We focused on the fiavor 
compositions of SUSY particle mass eigenstates, which refiect the relative size of fiavor-violating 
soft SUSY breaking terms. We enumerated three modes of collider experiments where one can 
measure the compositions by observing rare decays of SUSY particles. Predictions on their 
branching ratios were made based on the bulk matter RS model, and were compared with 
those of the minimal fiavor violation scenario. These predictions will be confirmed or rejected 
by a future lepton collider whose center-of-mass energy is tuned appropriately. 

A lesson of this study is that if new physics at TeV scale contains a fiavor-violating sector 
other than the Yukawa couplings, it is possible to observe signatures of models that explain 
the Yukawa coupling hierarchy through the fiavor structure of the new sector. In the case of 
this paper, MSSM contains gravity-mediation-origined soft mass terms, which provide a new 
source of fiavor violation. Gravity mediation and the Yukawa couplings are independent in 
the original MSSM, but have a correlation if the bulk matter RS model is the origin of the 
Yukawa coupling hierarchy. Hence we can predict the fiavor structure of gravity mediation 
contributions (up to their orders of magnitudes) from the data on SM, and eventually confirm 
or reject the bulk matter RS model through a detailed study on SUSY matter particles. This 
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study can be extended to any new physics scenario at TeV scale as long as it couples to matter 
fields and may violate flavor. In any case, SU(2) singlet muon and charm and their new physics 
partners play a pivotal role; SU(2) singlets receive less flavor-violating quantum corrections 
from the SM Yukawa couplings, and thus new flavor- violating terms are easy to extract. Since 
the 1st and 2nd generation particles only have small Yukawa couplings, we expect that their 
new SU(2) singlet partners almost do not mix with SU(2) doublets. Muon has much larger 
Yukawa coupling than electron and is much more sensitive to the origin of the Yukawa coupling 
hierarchy. The partner of charm may have a large flavor-violating mixing with top, the only 
quark whose flavor can be identified with virtually no misidentification rate. 
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